Colour field fiux tubes and Casimir scaling for various SU(3) representations 



N. Cardoso, M. Cardoso, P. Bicudo* 

CFTP, Departamento de Fisica, Instituto Superior Tecnico (Universidade Tecnica de Lisboa), Av. Rovisco Pais, 1049-001 

Lisboa, Portugal 



Abstract 

We investigate the QCD flux tubes linking static colour charges of different SU(3) representations, 
relevant to the understanding of confinement and of the Lund strings of Heavy Ion CoUisions. The colour field 
densities, the Casimir scaling factors and the widths for the flux tubes of the flrst flve different representations 
are computed in quenched SU(3) lattice QCD. This study is relevant to understand the mechanisms of 
conflnement and also the flux tubes utilized in the Lund Model and in other models of Heavy Ion Collisions. 
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1. Introduction 

We investigate the QCD flux tubes linking static 
colour charges of difl'erent SU(3) representations. 

Particle physics has been researching the build- 
ing blocks of nature, now compiled in a single Stan- 
dard Model, including the three fundamental quan- 
tum interactions. These - the strong, electromag- 
netic and weak interactions - together with gravi- 
tation (so far we don't know exactly how to quan- 
tize it) can all be described by simple mathematical 
equations, using the Lagrangian technique. Both 
the Standard Model and the gravitation Lagrangians 
have quite precise analytical solutions in simple ap- 
plications, say two-body problems or perturbative 
particle decays, except for the strong interaction. 
Because the internal symmetry group of the strong 
interaction is the non-commutative SU(3) , the strong 
interaction Lagrangian is too non-linear to have any 
analytical solution. Nevertheless the strong interac- 
tion - also denominated Quantum Chromodynam- 
ics (QCD) - Lagrangian may have numerical solu- 
tions when we utilize the Lattice QCD techniques 
devised by Wilson [1]. Lattice QCD utilizes the 
Dirac-Feynman Path integral, the Wick rotation, 
and the Boltzmann statistical mechanics to transfer 
the QCD problems into an average over a succession 
of configurations, similar to parallel but extremely 
small classical universes. Presently we are able to 
simulate lattices with sizes of the order of a few 
Fermi or femtometers. 

Confinement is one of the phenomena of QCD so 
far preventing its analytical solution. When a pair 
of QCD charge and anti-charge are slowly pulled 
apart, a flux tube of flelds develops between the 
pair. Flux tubes, sometimes approximated by a 
thin string for modelling, are extended and non- 
linear objects, and they have been observed in Lat- 
tice QCD, both with SU(2) and SU(3) gauge groups 
[2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13]. Although 
these flux tubes are not directly observed experi- 
mentally, since we don't know how to create static 
enough flux tubes in the Laboratory, their evidence 
is found in the spectrum of mesons composed of a 
quark, an antiquark, and the flux tube these QCD 
charges produce. 

A new type of QCD flux tubes, different from 
the mesonic flux tubes, may also be indirectly ob- 
served in the Heavy Ion Collisions presently pro- 
duced at the LHC and at the BNL. Strings or flux 
tubes are the building blocks of different models 
Heavy Ion Collision models. Heavy Ions, i. e. nu- 
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Table 1: Group theoretical factors for SU(3) as 
in [14]. D is the dimension of the representa- 
tion, (p, q) are the respective weight factors and 
z = exp(27r/3). We compare our colour field ratios 
with the casimir scaling factor djj = C'u / C3 , ratio 
of the quadratic Casimir charge of representation 
D and the one of the fundamental representation 3. 

clei composed of many particles, are collided at the 
highest possible energies, to approach in the Lab- 
oratory the extremely high temperatures occurring 
shortly after the Big Bang epoch in cosmological 
models. In the Lund hadronization model, ]15, 16, 
17] the large number of particles produced in Heavy 
Ion Collisions are reproduced with string fragmen- 
tation. Flux tubes also constitute the initial con- 
dition of recent hydrodynamic models describing 
the elliptic flow in Heavy Ion collisions [18]. Sur- 
prisingly, the flux tubes utilized to simulate Heavy 
Ion Collisions may have a width as thin as 0.2 fm, 
and a string tension almost one order of magnitude 
larger than the fundamental string tension linking a 
mesonic quark-antiquark pair. Moreover, the Colour 
Glass Condensate [19, 20[, produced at the onset 
of the Heavy Ion Collisions, is saturated by glu- 
ons. The quark degrees of freedom can be neglected 
at the onset of the quark gluon plasma, it is suffi- 
cient to study quenched SU(3) to address the Heavy 
Ion Collision flux tubes. It has also been proposed 
that the gluonic partons saturating the Heavy Ion 
Collisions produce in the Colour Glass Condensate 
perturbative flux tubes with an original width, or 
transverse size, of the order of l/Qs, with Q ~ 1 ^■ 
2 GeV [21, 22, 23[, flux tubes persisting during the 
evolution of the Quark Gluon Plasma. Thus it is 
relevant to study the conflning and non-perturbative 
flux tubes of higher SU(3) representations to un- 
derstand if they are compatible with the flux tubes 
produced at the Heavy Ion Collisions. 

The static potentials attracting a static colour 
charges and a colour anti-charge of different SU(3) 
representations have been already studied in quenched 
Lattice QCD. In the study of potentials, evidence 
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Figure 1: Lagrangian density created by static SU(3) charges of the representations, respectively from left 
to right: 3, 8, 6a, 15a, 10. The results are presented in lattice spacing units (colour online). 



was found for Casimir scaling in SU(2), first found 
by Anibjorn et al. [24, 25]. Casimir scaling occurs 
when the string a of the linear component of the 
static potential measured at intermediate distances 
is proportional to the Casimir invariant tr{A° • A"} 
was actually identified as and important topic by 
Del Debio et al. [26, 27[. In what concerns SU(3), 
Deldar [28, 29, 30[ and Bah [14[ have also found 
a good agreement with Casimir scaling. Note that 
the Casimir factors are computed analytically in 
group theory, see Table 1 for details. Bali, within 
an accuracy of 5 percent, found no violation to the 
Casimir scaling. To explain Casimir Scaling, Semay 
proposed that the cross section of the string is the 
same for all representations [31[. 

Here we observe the colour fiux tubes produced 
by a charge and an anti-charge of higher SU(3) rep- 
resentations. We extend the techniques developed 
by the previous authors both to study fiux tubes 
and static potentials of higher SU(3) representa- 
tions. Previously only the lowest 3 representation, 
and more recently the 8 representation were studied 
[32, 33[. 

In this paper we present the results for the colour 
fields, the energy and Lagrangian densities as well 
as the respective Casimir scaling factors, using the 
energy density along the fiux tube, for five differ- 
ent representations, computed in quenched SU(3) 



lattice QCD. In section II, we introduce the lattice 
QCD formulation. We briefly review the Wilson 
loops for this system, which was used in [14[, and 
show how we compute the colour fields as well as 
the Lagrangian and energy density distribution. In 
section III the numerical results are shown, includ- 
ing the chromoelectric field, chromomagnetic field, 
Lagrangian and energy densities, the Casimir scal- 
ing factors and the fiux tube widths for the various 
SU(3) representations. Finally, we present our con- 
clusion in section IV. 

2. The Wilson Loops and Colour Fields 

To impose the position of the static charge and 
anti-charge in our Lattice, we utilize a Wilson loop. 
The Wilson loop W£)(ri,r2,T) creates a charge at 
position Ti and an anti-charge at Vi separated by 
a spatial distance i? = |r2 — ri|, which propagate 
for an Euclidian time T maintaining their spatial 
positions ri and r2unchanged. To preserve gauge 
invariance, the wilson loop consists of a simple rect- 
angular loop of size RxT composed of the trace of 
the product of links Ud of the respective represen- 
tation D of SU(3). 

Utilizing the algebra of the Gell-Mann gener- 
ators of SU(3), the Wilson loops W^D(ri,r2,T) for 
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Figure 2: Chromoelectric field, chromomagnetic 
field, Lagrangian and energy densities produced by 
the various SU(3) charges in the y axis {x = z = 0). 
The charges are located at y = ±4. 



Figure 3: Chromoelectric field, chromomagnetic 
field, Lagrangian and energy densities produced by 
the various SU(3) charges along a radial distance in 
the metriatrix plane y = 0. 



the various SU(3) representations [14], can be trans- 
lated into products of the fundamental representa- 
tion D = 3 Wilson loops W3 , and fundamental links 



1/3 = 11. U denotes a group element in the funda- 
mental representation of SU(3). The traces of Ud 
in various representations, can easily be expressed 
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in terms of traces of powers of U as in Bali et al. 
[14], 



= Tr [/ , 
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The chromoelectric and clironioniagnetic fields 
on the lattice are given by the Wilson loop and 
plaquette expectation values, 
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where the jk indices of the plaquette complement 
the index i of the magnetic field, and where the 
plaquette at position r = {x,y,z) is computed at 
t = T/2, 

P^u (r) - 1-^Re Tr [U^ir)U,ir + fi)Ul{r + ly)Ut(r)] 

(3) 

Importantly, in SU(3) lattice QCD, since the Gel- 
Mann matrices have zero trace, the plaquette only 
produces the square of the components of the fields. 
This is consistent with the fact that the colour field 
components are not gauge invariant (unlike the fields 
in electrodynamics). 

Nevertheless in SU(3) we can compute the La- 
grangian (£) and energy (H) densities, since they 
are computed from the the square of the field com- 
ponents, 

(H(r)) = i((E2(r)) + (B2(r))) , (4) 
(£(r)> = i((E2(r))-(B2(r))) . (5) 

To compute the static field expectation value, 
we plot the expectation value {Ef{r)') or (i3^(r)), 
as a function of the temporal extent T of the Wilson 
loop. At sufficiently large T, the groundstate cor- 
responding to the studied quantum numbers dom- 
inates, and the expectation value tends to a hori- 
zontal plateau. 

In order to improve the signal to noise ratio of 
the Wilson loop, we use the APE smearing defined 



by 

+w J2 Uu {s) {s + v) Ul (s + m) ()6) 

with w ~ 0.2 and iterate this procedure 50 times 
in the spatial direction [32]. To achieve better ac- 
curacy in the flux tube, we apply three levels of 
hypercubic blocking (HYP) in the time direction, 
[34], with 



ai — 0.75, a2 — 0.6, — 0.3 



(7) 



Note that these two procedures are only applied to 
the Wilson Loop, not to the plaquette. To compute 
the fields, we fit the horizontal plateaux obtained 
for each point (x, y, z) determined by the plaquette 
position, but due to the azimuthal symmetry we 
only consider two coordinates for simplicity. For 
the positions x (distance to the charge axis) and 
y (distance along the charge axis) considered, ex- 
cept at the charge location, we find in the range 
of T G [3, 12] in lattice units, horizontal plateaux 
with a /dof e [0.3,2.0]. We find these hori- 
zontal plateaux for representations up to Z? = 10. 
For higher representations then the statistical noise 
increases significantly already at T below 12, and 
thus we decide not to further study representations 
higher than D = 10. We finally compute the error 
bars of the fields with the jackknife method. 

3. Results 

Here we present the results of our simulations 
with 381 quenched configurations in a 32^ X 64 lat- 
tice at P = 6.2. We generate our configurations in 
NVIDIA CPUs of the FERMI series (480, 580 and 
Tesla 2070) with a SU(3) CUBA code upgraded 
from our SU(2) combination of Cabibbo-Marinari 
pseudoheatbath and over-relaxation algorithm [35, 
36]. Our SU(3) updates involve three SU(2) sub- 
groups, we work with 9 complex numbers, and we 
reunitarize the matrix. We have two options to save 
the configurations, either in a structure of arrays 
where each array lists a given complex component 
for all the lattice sites, or in an array of structures 
where each structure is a SU(3) matrix. 

The results are presented in lattice spacing units 
of a, with a = 0.07261(85) fm or a"! = 2718 ± 32 
MeV. We locate our colour charge and anti-charge 
at ri = (0, -4, 0) and ra = (0, 4, 0). 
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Figure 4: Casimir scaling of the energy density 
for the various SU(3) representations in the charge 
axis. 
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Figure 5: Casimir scaling of the energy density for 
the various SU(3) representations in the mediator 
plane. 



3.1. Colour Fields 

In Figure 1 we depict the 3D plots of the flux 
tubes produced by the various SU(3) charges. 

To arrive at a more quantitative picture of the 
flux tubes, we cut the flux tubes in their symme- 
try axis. We depict in Figure 2 the proflles of the 
different components of the flux tube, in a plane 
parallel to the quark-antiquark axis. It is is well 
known that the APE smearing and HYP blocking 
affect the fields in the charge neighbourhood [37], 
decreasing the respective peaks and shifting them 
slightly from the ri and r2, but do not affect the aim 
of our study which is the study of the flux tubes. 

We also show in Figure 3 the profiles of the dif- 
ferent components of the fiux tube, in the mediator 
plane, perpendicular to the quark-antiquark axis. 

The profiles of the different representations are 
similar in shape, only changing in magnitude, as 
suggested by Semay to understand Casimir Scaling 
[31]. 

3.2. Casimir Scaling 

We further compare the different SU(3) flux tubes 
with the Casimir scaling. 

In the Figs. 4 and 5, we show the results for 
the ratio between the energy density of every rep- 
resentation over the energy density of the funda- 
mental triplet representation. We compare these 
ratios with the ratio of the Casimir scaling hypoth- 
esis, is correct. Indeed, at least for the first five 
lower representations of SU(3), and in the centre of 



the flux tube where most of the colour fleld inten- 
sity is concentrated, our results are in agreement 
with a flux tube Casimir scaling, within error bars. 

3.3. Flux tube width 

Moreover, we determine the width of the vari- 
ous SU(3) flux tubes, computing the Radius Mean 
Square in the mediator plane ?; = 0, utilizing the 
Lagrangian density as a probability density. 



RMS 



(8) 



evaluated only in the positive x direction. In Ta- 
ble 2 we show the respective RMS obtained for 
the various SU(3) representations. We find a simi- 
lar Radius Mean Square for the various SU(3) flux 
tubes we study, of the order of 0.35 fm. This is 
consistent with the conjecture of Semay ]31] for a 
constant cross section of the various flux tubes. 

4. Conclusions 

We study in quenched Lattice QCD the flux 
tubes produced by static SU(3) colour charges. This 
is relevant to understand the mechanisms for the 
confinement of colour SU(3) charges, and also to 
find a solid foundation for the very strong and thin 
flux tubes utilized in models of Heavy Ion Colli- 
sions. 

For the flrst flve representations 3, 8, 6, 15a and 
10, up to Casimir scaling factors of 4.5, our Wilson 
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_D BMS lliii] 

~3 0.351 ±0.023 

8 0.343 ± 0.028 

6 0.392 ± 0.029 

15« 0.312 ±0.054 

10 0.395 ±0.123 



Table 2: The different Radius Mean Squares ob- 
tained for each of the SU(3) representations we 
study here. 

loop plateaux have acceptable /dof S [0.3,2.0] 
and we depict the different colour field components. 

Our computations show that the energy density 
of the flux tubes produced by the different static 
SU(3) sources compatible with Casimir scaling. The 
energy density, within error bars, is proportional to 
the eigenvalue of the quadratic Casimir operator of 
that representation. 

All the flux tubes of the different representations 
show a similar radius mean squares circa 0.35 fm. 
This is twice as large, but still of the same order, of 
the width of the thin flux proposed to exist at the 
onset of Heavy Ion Collisions. It remains possible 
that the so-called perturbative flux tubes of Heavy 
Ion Collisions evolve to the non-perturbative flux 
tubes observed here. 

Acknowledgments 

We thank Edmond lancu, Gunnar Bali, Jeff Green- 
site, Magdalena Malek and Orlando Oliveira for 
useful discussions. 

This work was partly funded by the FCT pluri- 
annual contract, PTDC/FIS/100968/2008, and the 
FCT annual contracts CERN/FP/109327/2009 and 
CERN/FP/116383/2010. Nuno Cardoso is supported 
by the FCT fellowship SFRH/BD/44416/2008. 

References 

[1] K. G. Wilson, Confinement of Quarks. Phys.Rcv. DIO 
(1974) 2445-2459. 

[2] M. Pukugita, T. Niuya, DISTRIBUTION OF CHRO- 
MOELECTRIC FLUX IN SU(2) LATTICE GAUGE 
THEORY, Phys.Lett. B132 (1983) 374. 

[3] J. W. Flower, S. W. Otto, THE FIELD DISTRI- 
BUTION IN SU(3) LATTICE GAUGE THEORY, 
Phys.Lett. B160 (1985) 128. 

[4] R. Sommcr, CHROMOFLUX DISTRIBUTION IN 
LATTICE QCD, NucI.Phys. B291 (1987) 673-691. 



[5] D. Caldi, T. Sterling, GRAPHIC EVIDENCE FOR 
NONABELIAN FLUX TUBES, Phys.Rev.Lett. 60 
(1988) 2454. 

[6] A. Di Giacomo, M. Maggiore, S. Olejnik, EVIDENCE 
FOR FLUX TUBES FROM COOLED QCD CONFIG- 
URATIONS, Phys.Lett. B236 (1990) 199. 

[7] A. Di Giacomo, M. Maggiore, S. Olejnik, CONFINE- 
MENT AND CHROMOELECTRIC FLUX TUBES IN 
LATTICE QCD, NucI.Phys. B347 (1990) 441-460. 

[8] M. S. Cardaci, P. Cea, L. Cosmai, R. Falcone, A. Papa, 
Chromoelectric flux tubes in QCD, Phys.Rev. D83 
(2011) 014502. 

[9] C. Michael, Adjoint Sources in Lattice Gauge Theory, 
NucI.Phys. B259 (1985) 58. 
[10] H. Markum, M. Faber, QUARK MULTIPLETS AND 
THEIR INTERACTIONS, Phys.Lett. B200 (1988) 
343-347. 

[11] N. Campbell, I. Jorysz, C. Michael, THE ADJOINT 
SOURCE POTENTIAL IN SU(3) LATTICE GAUGE 
THEORY, Phys.Lett. B167 (1986) 91. 

Jl2] G. Bali, K. Schilling, C. Schlichter, Observing long color 
flux tubes in SU(2) lattice gauge theory, Phys.Rev. D51 
(1995) 5165-5198. 

Jl3] F. Bissey, A. I. Signal, D. B. Leinweber, Compari- 
son of gluon flux-tube distributions for quark- diquark 
and quark-antiquark hadrons, Phys. Rev. D80 (2009) 
114506. 

[14] G. S. Bali, Casimir scaling of SU(3) static potentials, 

Phys. Rev. D62 (2000) 114503. 
[15] B. Andersson, G. Gustafson, B. Soderberg, A General 

Model for ,let Fragmentation, Z.Phys. C20 (1983) 317. 
]16] B. Andersson, G. Gustafson, G. Ingelman, T. Sjos- 

trand, Parton Fragmentation and String Dynamics, 

Phys.Rept. 97 (1983) 31-145. 
]17] B. Andersson, The Lund model, 

Camb.Monogr. Part. Phys. Nucl. Phys. Cosmol. 7 (1997) 

1-471. 

]18] P. Mota, T. Kodama, T. Koide, Multi-flux tube ini- 
tial condition and event-by-event hydrodynamics, Nucl. 
Phys. A862-863 (2011) 188 191. 

]19] E. lancu, A. Lcoiiidov, L. D. McLerrari, Nonlinear gluon 
evolution in the color glass condensate. I, Nucl. Phys. 
A692 (2001) 583-645. 

]20] E. lancu, A. Leonidov, L. D. McLerran, The renormal- 
ization group equation for the color glass condensate, 
Phys. Lett. B510 (2001) 133-144. 

]21] A. Dumitru, F. Gelis, L. McLerran, R. Venugopalan, 
Glasma flux tubes and the near side ridge phenomenon 
at RHIC, Nucl. Phys. A810 (2008) 91-108. 

]22] K. Dusling, F. Gelis, T. Lappi, R. Venugopalan, Long 
range two-particle rapidity correlations in A+A colli- 
sions from high energy QCD evolution, Nucl. Phys. 
A836 (2010) 159-182. 

]23] T. Lappi, L. McLerran, Long range rapidity correlations 
as seen in the STAR, experiment, Nucl. Phys. A832 
(2010) 330-345. 

]24] J. Ambjorn, P. Olesen, C. Peterson, Stochastic Confine- 
ment and Dimensional Reduction. 1. Four-Dimensional 
SU(2) Lattice Gauge Theory, NucI.Phys. B240 (1984) 
189. 

]25] J. Ambjorn, P. Olesen, C. Peterson, STOCHAS- 
TIC CONFINEMENT AND DIMENSIONAL REDUC- 
TION. 2. THREE-DIMENSIONAL SU(2) LATTICE 
GAUGE THEORY, NucI.Phys. B240 (1984) 533. 

]26] L. Del Debbio, M. Faber, J. Greensite, S. Olejnik, 



7 



Casimir scaling versus Abelian dominance in QCD 

string formation, Phys.Rev. D53 (1996) 5891-5897. 
[27] L. Del Debbio, M. Faber, J. Greensite, S. Olejnik, Some 

cautionary remarks on Abelian projection and Abelian 

dominance, Nucl. Phys. Proc. Suppl. 53 (1997) 141-147. 
[28] S. Deldar, The string tensions of SU(3) representations, 

Nucl. Phys. Proc. Suppl. 73 (1999) 587-589. 
]29] S. Deldar, Static SU(3) potentials for sources in various 

representations, Phys. Rev. D62 (2000) 034509. 
]30] S. Deldar, A new lattice measurement for potentials 

between static SU(3) sources, Eur. Phys. J. C47 (2006) 

163-169. 

]31] C. Semay, About the Casimir scaling hypothesis, Eur. 
Phys. J. A22 (2004) 353-354. 

]32] M. Cardoso, N. Cardoso, P. Bicudo, Lattice QCD com- 
putation of the colour fields for the static hybrid quark- 
gluon-antiquark system, and microscopic study of the 
Casimir scaUng, Phys. Rev. D81 (2010) 034504. 

]33] N. Cardoso, M. Cardoso, P. Bicudo, Gauge invariant 
SU(3) lattice computation of the dual gluon mass and of 
the dual Ginzburg-Landau parameters A and t; in QCD 

(2010) . 

]34] A. Hasenfratz, F. Knechtli, Flavor symmetry and the 
static potential with hypercubic blocking, Phys. Rev. 
D 64-3 (2001) 034504. 

]35] N. Cardoso, P. Bicudo, SU(2) Lattice Gauge Theory 
Simulations on Fermi GPUs, J. Comput. Phys. 230 

(2011) 3998-4010. 

]36] 2011. The CUDA codes are available at Portuguese 
Lattice QCD collaboration, http://nemea.ist.utl. 
pt/''ptqcd. 

]37] F. Okiharu, R. M. Woloshyn, An alternate smearing 
method for Wilson loops in lattice QCD, Eur. Phys. J. 
C35 (2004) 537-542. 



8 



